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ABSTRACT 

The R CrA region was observed in the 3.5 and 6.2 cm continuum with high angular resolutions (0.6-1.7 arcseconds) 
using the Very Large Array. Archival data sets were also analyzed for comparison, which provided angular resolutions 
up to 0.3 arcseconds. A cluster of young stellar objects was detected, and a rich array of star forming activities was 
revealed. IRS 7A showed an enhanced outflow activity recently. The main peak of IRS 7A positionally coincides with 
an X-ray source, which suggests that the X-ray emission is directly related to the central protostar The Class source 
SMA 2 is associated with a double radio source, B 9a and 9b, and seems to be driving two outflows. The B 9 complex 
is probably a multiple-protostar system. Both B 9a and 9b are nonthermal radio sources with negative spectral indices. 
IRS 7B is a compact radio source surrounded by an extended structure. The compact source corresponds to the Class O/I 
source SMA 1, and it is also closely associated with an X-ray source, suggesting that magnetic activities start early in 
the protostellar stage of evolution. The extended structure of IRS 7B may be a bipolar outflow. IRS 5 was resolved into 
two sources with a separation of 0.9 arcseconds. Both IRS 5a and 5b display radio flares and X-ray emission, suggesting 
that energetic magnetic processes are active in both members. The month-scale active phase of IRS 5b implies that the 
flare activity must involve large-scale magnetic fields. During the strong flare event of IRS 5b in 1998, IRS 5a also 
showed an enhanced level of radio emission. This concurrent activity suggests that IRS 5 may be an interacting young 
binary system, but the interaction mechanism is unknown. Alternatively, what was seen in the radio images could be a 
circumbinary halo. The variable radio source B 5 was found to be a nonthermal source at times, and the size of the 6.2 
cm source is about 1 arcsecond, suggesting that B 5 is a Galactic object. A radio outburst of IRS 6 was detected once, 
and the radio/X-ray source was identified as IRS 6a. The other member of the IRS 6 system, IRS 6b, was undetected in 
X-rays, suggesting that only IRS 6a has detectable magnetic activities. Properties of other radio sources, IRS 1, IRS 2, 
and R CrA, are discussed, and the radio detections of T CrA and WMB 55 are reported. The proper motion of R CrA 
was marginally detected. Also presented is the classification of infrared sources in the R CrA region based on an infrared 
color-color diagram. 



Subject headings: binaries: general — 
ISM — stars: flare 

1. INTRODUCTION 



ISM: individual (R Coronae Australis) 
— stars: formation 



ISM: structure — radio continuum: 



The Corona Australis molecular cloud is one of the nearest 
star-forming regions, and the R CrA region is the most active 
site of star formation in this cloud (Harju et al. 1993; Chini et 
al. 2003). The R CrA cloud core contains a cluster of young 
stellar objects known as the "Coronet" cluster (Taylor & Storey 
1984; Wilking et al. 1992, 1997). Deeply embedded protostars 
and preprotostellar cores were detected in the (sub)millimeter 
continuum (Henning et al. 1994; Choi & Tatematsu 2004; Nut- 
ter et al. 2005; Groppi et al. 2004, 2007). Star forming activities, 
such as infall, accretion, molecular outflows, and Herbig-Haro 
objects, were revealed by millimeter spectroscopic observations 
and optical imaging (Strom et al. 1974; Hartigan & Graham 
1987; Anderson et al. 1997; Groppi et al. 2004, 2007; Nisini 
et al. 2005; Wang et al. 2004). Centimeter continuum imaging 
with interferometers has been especially effective in probing the 
nature of the deeply embedded objects and determining their 
precise positions (Brown 1987; Suters et al. 1996; Feigelson 
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et al. 1998; Forbrich et al. 2006; Miettinen et al. 2008). High- 
energy processes around some of the young stellar objects were 
discovered in X-ray observations (Koyama et al. 1996; Ham- 
aguchi et al. 2005, 2008; Forbrich et al. 2006, 2007; Forbrich 
& Preibisch 2007). The distance to the CrA molecular cloud is 
-170 pc (Knude & H0g 1998). 

Among the near-IR sources discovered by Taylor & Storey 
(1984) in the R CrA cloud core, IRS 7 is the most deeply em- 
bedded and probably the youngest object. In the centimeter con- 
tinuum imaging by Brown (1987), IRS 7 was resolved into three 
radio sources: IRS 7A, IRS 7B, and B 9. Single-dish observa- 
tions of submillimeter continuum emission suggested that they 
are deeply embedded objects (Nutter et al. 2005). Several ra- 
dio continuum sources were detected around IRS 7A (Choi & 
Tatematsu 2004; Forbrich et al. 2006), suggesting that this small 
region contains a multiple-protostar system. The IRS 7A region 
is interesting because it appears to contain objects with inho- 
mogeneous spectral classes or in different stages of evolution 
(Choi & Tatematsu 2004). While IRS 7A appears to be a Class 
I protostar as it was detected in near-IR, the other members of 
the multiple system (CT 3 and B 9) were suggested to be prob- 
able Class protostars. Recent interferometric observations in 
the millimeter continuum confirmed that B 9 (SMA 2) is a Class 
source while IRS 7A and CT 3 were undetected (Groppi et al. 
2007). 

IRS 7B attracted attention because it is a rare example of 
Class protostars showing X-ray emission (Feigelson & Mont- 
merle 1999; Hamaguchi et al. 2005). Recent interferometric 
observations in the millimeter continuum suggested that IRS 
7B (SMA 1) may be a Class O/I transitional object (Groppi et 
al. 2007). Hamaguchi et al. (2005) demonstrated that the X- 
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Table 1 

Parameters of the VLA Observing Runs 



Synthesized Beam 











Frequency 


Natural weight 


Uniform weight 


Array 


Track 


Date 




(GHz) 


Size 


P.A. 


Size 


P.A. 


Configuration 


Tr 1 


1996. 


12. 


29. 


4.9 


1.21 


X 


0.50 


-6.5 


0.86 


X 


0.35 


—7.1 


A 










8.5 


0.70 


X 


0.31 


-2.5 


0.51 


X 


0.20 


0.1 


A 


Tr 2 


1997. 


01. 


19. 


8.5 


0.82 


X 


0.59 


10.6 


0.57 


X 


0.20 


0.8 


BnA 


Tr 3 


1997. 


01. 


20. 


8.5 


1.04 


X 


0.89 


— 10.4 


0.84 


X 


0.22 


10.6 


BnA 


Tr 4 


1998. 


06. 


27. 


8.5 


1.13 


X 


0.69 


23.8 


0.82 


X 


0.52 


41.0 


BnA 


Tr 5 


1998. 


07. 


19. 


8.5 


2.61 


X 


0.78 


-1.9 


2.38 


X 


0.54 


—0.8 


B 


Tr 6 


1998. 


09. 


07. 


8.5 


2.61 


X 


0.78 


5.4 


2.37 


X 


0.54 


4.7 


B 


Tr 7 


1998. 


09. 


19. 


8.5 


2.54 


X 


0.80 


0.0 


2.30 


X 


0.53 


0.4 


B 


Tr 8 


1998. 


09. 


27. 


8.5 


2.71 


X 


0.79 


11.6 


2.43 


X 


0.52 


11.8 


B 


Tr 9 


1998. 


10. 


02. 


8.5 


2.58 


X 


0.79 


-5.4 


2.25 


X 


0.54 


-4.6 


B 


Tr 10 


1998. 


10. 


06. 


8.5 


2.65 


X 


0.79 


-6.6 


2.43 


X 


0.54 


-5.2 


B 


Tr 11 


1998. 


10. 


10. 


8.5 


2.54 


X 


0.79 


0.1 


2.30 


X 


0.53 


0.4 


B 


Tr 12 


1998. 


10. 


13. 


8.5 


2.60 


X 


0.80 


-5.6 


2.25 


X 


0.54 


-4.6 


B 


Tr 13 


2005. 


02. 


03. 


8.5 


1.01 


X 


0.74 


23.4 


0.73 


X 


0.54 


33.9 


BnA 


Tr 14 


2005. 


02. 


12. 


4.9 


2.10 


X 


1.35 


9.9 


1.45 


X 


1.02 


14.7 


BnA 



Note. — Units of beam size and P.A. are arcseconds and degrees, respectively. Tracks Tr 2-3 correspond to the observing 
runs presented by Feigelson et al. 1998. Tracks Tr 4—12 correspond to epochs Rl-9 in Table 1 of Forbrich et al. 2006. 



ray emission from IRS 7B comes from two different compo- 
nents: one hot and variable, and the other cool and constant. The 
cool component was interpreted as originating from collisional 
shock of the protostellar jet with circumstellar gas, and such X- 
ray emission has been found around several Class O/I protostars, 
for example, HH 2, L1551 IRS 5, and OMC 3 MMS 2 (Pravdo 
et al. 2001; Favata et al. 2002; Bally et al. 2003; Tsujimoto et al. 
2004). The hot variable X-ray emission had not been detected 
from other Class protostars. This component would be related 
to the reconnection of magnetic fields around the protostar dur- 
ing mass accretion (e.g., Montmerle et al. 2000). Spatially sep- 
arating the jet collision region from the central protostar with a 
high resolution imaging is therefore important in understanding 
the structure of the protostellar system. 

IRS 5 is an interesting young binary system. It was spatially 
resolved in the near-IR and X-ray images (Nisini et al. 2005; 
Hamaguchi et al. 2008). Both members of the system show qui- 
escent X-ray emission, but flare activities of them were remark- 
ably different: only IRS 5 a frequently showed solar type X-ray 
flares (Hamaguchi et al. 2008). IRS 5 also displays radio out- 
bursts and emits circularly polarized radio emission (Suters et 
al. 1996; Feigelson et al. 1998; Forbrich et al. 2006; Miettinen 
et al. 2008). However, the radio properties of each binary mem- 
ber need to be characterized separately to understand the radio 
activities. 

The R CrA region also contains other radio sources display- 
ing interesting star-formation activities of their own. IRS 1 and 
IRS 2 are Class I protostars and variable X-ray sources (For- 
brich et al. 2006). IRS 1 is the driving source of the HH 100 
outflow (Strom et al. 1974). WMB 55 is a Class I source with- 
out an X-ray detection (Wilking et al. 1997; Nutter et al. 2005). 
IRS 6 is probably an X-ray emitting T Tauri star driving a giant 
Herbig-Haro flow (Wang et al. 2004; Forbrich et al. 2006). R 
CrA and T CrA are Herbig Ae/Be stars, and each of them may 
be a binary system (Takami et al. 2003). B 5 is a radio source of 
unknown nature, and its flux varies rapidly (Brown 1987; Suters 



et al. 1996; Feigelson et al. 1998). 

To obtain high-quality images of the young stellar objects 
despite the low declination, we observed the R CrA region with 
the Very Large Array (VLA) using a fast-switching phase cali- 
bration technique. In this paper, we present our centimeter con- 
tinuum observations of the R CrA region. We describe our radio 
continuum observations in § 2 and archival data in § 3. In § 4 
we briefly report the results of the continuum imaging. In § 5 
we discuss the star forming activities of the radio sources in the 
R CrA region. A summary is given in § 6. 

2. OBSERVATIONS 

The R CrA IRS 7 region was observed using VLA of the 
National Radio Astronomy Observatory (NRAO) in the stan- 
dard X-band continuum mode (8.5 GHz or A = 3.5 cm) and 
in the standard C-band continuum mode (4.9 GHz or A = 6.2 
cm). Twenty-three antennas were used in the BnA-array config- 
uration on 2005 February 3 for the 3.5 cm observations and on 
2005 February 12 for the 6.2 cm observations (Tr 13/14 in Ta- 
ble 1). The phase-tracking center was a2ooo = 19''01'"55.33s, 
(^2000 = -36°57'21.66". The phase was determined by observ- 
ing the nearby quasar 1924-292 (QSO B1921-293). The target- 
calibrator cycle was 4 and 5 minutes for the 3.5 and 6.2 cm 
observations, respectively. In the 3.5 cm observations, the flux 
was calibrated by setting the flux density of 1924-292 to 13.7 
Jy, which is the flux density measured within 16 days of our ob- 
servations (VLA Calibrator Flux Density Databaso- To check 
the flux scale, the quasar 1331H-305 (4C 30.26) was observed in 
the same track. Comparison of the amplitude gave a flux den- 
sity of 5.30 Jy for 1331H-305, which agrees with the value in the 
VLA Calibrator Flux Density Database to within 3%. In the 6.2 
cm observations, the flux was calibrated by setting the flux den- 
sity of 1331H-305 to 8.25 Jy, which is the flux density measured 
within 2 days of our observations according to the VLA Calibra- 

'See |http://aips2.nrao.edu/vla/calflux.htmr| 
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tor Flux Density Database. Comparison of the amplitude gave 
a flux density of 13.52 Jy for 1924-292. Because of the short 
cycle time, the quality of phase calibration was good. We tried 
a self-calibration, but it did not make a noticeable difference 
in the quality of resulting images. Therefore, the Tr 13/14 im- 
ages presented in this paper were made without self-calibration. 
Maps were made using a CLEAN algorithm. 

3. ARCHIVAL DATA 

In addition to our data, we also analyzed several VLA 
data sets retrieved from the NRAO Data Archive System. We 
searched for 3.5/6.2 cm continuum observations of the R CrA 
region made in the A-, BnA-, or B-array configuration and 
found four observing projects in recent years. Including our ob- 
servations, we present the results of 14 observing runs. A sum- 
mary of the observing runs is given in Table 1. The archival 
VLA data were obtained from observations with a typical 
target-calibrator cycle time of 10-15 minutes. To improve the 
phase calibration we used the phase-only self-calibration tech- 
nique (see Forbrich et al. 2006). The Tr 1-12 images presented 
in this paper were made with the self-calibration applied. 

3.1. Project AB 807 

The NRAO observing project AB 807 includes a single ob- 
serving run (Tr 1). The data set contains data from the observa- 
tions of the R CrA IRS 7 region in the standard X-band contin- 
uum mode and in the standard C-band continuum mode in a sin- 
gle track. Twenty-five antennas were used in the A-array config- 
uration. The phase-tracking center was 01950 = 18'^58™33.0s, 
<^i950 = -37°01'41.0". The phase caUbrators were the quasars 
QSO B 1921-293 and QSO B 1933^00, and the flux caUbrator 
was the quasar 3C48. The bootstraped flux densities were 17.4 
Jy for QSO B 1921-293 and 0.96 Jy for QSO B 1933^00 in the 
X-band and 16.2 Jy for QSO B 1921-293 and 0.98 Jy for QSO 
B1933-400 in file C-band. 

3.2. Project AC 481 

The NRAO observing project AC 481 includes two observ- 
ing runs (Tr 2 and 3). Details of the observations and results 
were presented by Feigelson et al. (1998). The flux densities of 
sources in the images produced by us agree well with the values 
reported by Feigelson et al. (1998). However, the positions of 
compact sources have an offset compared with those found in 
the images of the other data sets. The amount of the position 
offset is 0. 1" to the east and 0.4" to the north. A probable cause 
of this offset is that the phase calibrator used in the project AC 
481 was 1744-312. In the other projects the phase calibrator 
was 1924-292, as in our observations. The position accuracy of 
1744-312 is worse than that of 1924-292: 1744-312 is a code- 
C calibrator while 1924-292 is a code- A calibrator. The angular 
distance from the target source is also larger: 17.5° for 1744- 
312 and 9.2° for 1924-292. Therefore, we conclude that the 
astrometry of the AC 48 1 images is less reliable than that of the 
other images. In all the figures presented in this paper, a position 
correction was applied to the AC 481 images. 

3.3. Projects AK 469 and AM 596 

The NRAO observing projects AK 469 and AM 596 include 
nine observing runs (Tr 4-12). Details of the observations and 
results were presented by Forbrich et al. (2006). The flux densi- 
ties and positions of compact sources in the images produced by 
us agree weU with the values reported by Forbrich et al. (2006). 



3.4. SpitzerData 

In addition to the radio data, we analyzed archival infrared 
data to classify the young stellar objects in the R CrA region. 
From the Spitzer data archive, GTO observations for the R CrA 
region (AOR keys: IRAC 3650816 and MIPS 3664640) were 
retrieved. The aperture photometry was performed with mosaic 
post-BCD images, which were processed by the SSC pipeline 
(Ver. S14 and S16 for IRAC and MIPS, respectively), at aU 
IRAC bands (3.6, 4.5, 5.8, and 9.0 fim) and MIPS 24 and 70 
fim using the SSC Mopex package. The short-exposure high- 
dynamic range observations were used for the IRAC photome- 
try. 

4. RESULTS 

Table 2 lists aU the 16 sources detected in the radio maps. 
Figure la shows the average map of the R CrA region in the 
3.5i cm continuum from all the data sets, and Figure lb shows 
the maps of the IRS 7 region from our observations. All the im- 
ages presented in this paper are corrected for the primary beam 
response, and the lowest contour level roughly corresponds to 
the 3a value, where a is the rms of noise. Objects associated 
with the centimeter continuum sources are marked in the fig- 
ures, and the size of markers roughly corresponds to the posi- 
tion uncertainty of the objects. Tables 2 and 3 list the parameters 
of detected sources, and Figure 2 shows the flux variability of 
bright sources. We briefly summarize the overall properties of 
detected sources in this section, and detailed discussion on in- 
dividual sources will be given in § 5. 

The IRS 7 region shows a complex of several compact 
sources and some extended emission structures. Since there 
are several pubUshed papers reporting the centimeter contin- 
uum sources in this region (many of them have their own list 
of sources), simply using source numbers can be confusing. In 
this paper, the source names listed in Table 3 are used. We tried 
to list the most commonly used names. In some cases, for ex- 
ample, a name may be a combination of an acronym (names of 
authors who identified the source first) and a sequence (source 
number) used by those authors. 

4.1. Source Structure 

IRS 7A is the strongest centimeter-continuum source in this 
region and is composed of a compact source and an extended 
emission structure. B 9 is associated with the 7 nmi source CT 
2 and the 1 . 1 mm source SMA 2 and is probably the most deeply 
embedded source in this region. B 9 was resolved into two com- 
pact sources. FPM 10 and FPM 13 are extended emission struc- 
tures, and their shapes and peak positions change depending on 
the uv coverage and the epoch of observations. 

IRS 7B is also composed of a compact source surrounded 
by an extended emission structure. The compact source is asso- 
ciated with the 1 . 1 mm source SMA 1 . The extended structure 
clearly has a bipolar morphology elongated in the northeast- 
southwest direction. The northeastern lobe corresponds to FPM 
15 (Forbrich et al. 2006) and shows a clumpy structure. 

IRS 5 was resolved into two compact sources, and they dis- 
played interesting flux variability. B 5 is slightly extended in the 
northeast-southwest direction, and the deconvolved size along 
the elongation is ^1". Most of the other detected sources in the 
R CrA field are point-like objects. 

4.2. Spectral Energy Distribution 

Table 3 lists a, the spectral index in the centimeter band. For 
most sources, the spectral index is consistent with the value ex- 
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Table 2 

Centimeter Continuum Sources in the R CrA Region 



Peak Position'' 

Source aj2ooo.o 5j2ooo.o Associated Objects'' 

1 19 01 41.58 -36 58 31.0 IRS 2, FCW 1, SMM 5 

2 19 01 43.27 -36 59 12.3 B5 

3 19 01 48.02 -36 57 22.5 IRS 5a, B 7, SMM 4, IRAC 10 

4 19 01 48.06 -36 57 21.8 IRS 5b, B 7, SMM 4 

5 19 01 50.46 -36 56 37.6 IRS 6, FPM 6, IRAC 8 

6 19 01 50.67 -36 58 09.7 IRS 1, B 8, SMM 3, IRAC 7 

7 19 01 53.68 -36 57 08.2 R CrA, FCW 5, IRAC 6 

8 19 01 55.04 -36 57 16.6 FPM 10 

9 19 01 55.26 -36 57 16.9 B 9, CT 2, SMM IC, SMA 2 

10 19 01 55.30 -36 57 16.7 B 9, CT 2, SMM IC, SMA 2 

11 19 01 55.32 -36 57 22.1 IRS 7A, CT 4, IRAC 5 

12 19 01 55.34 -36 57 12.9 FPM 13 

13 19 01 55.61 -36 57 43.4 SMM lA 

14 19 01 56.41 -36 57 28.1 IRS 7B, SMM IB, SMA 1, IRAC 4 

15 19 01 58.55 -36 57 08.5 WMB 55, SMM 2, IRAC 3 

16 19 01 58.78 -36 57 49.9 T CrA, IRAC 1 



Note. — Units of right ascension are hoiirs, minutes, and seconds, and units of 
declination are degrees, arcminutes, and arcseconds. 

* For most sources, the positions are from the natural-weight 3.5 cm continuum 
map of Tr 1. For source 1 (IRS 2), the position is from the 6.2 cm map of Tr 1. For 
source 5 (IRS 6), the position is from the 3.5 cm map of Tr 12. For sources 8, 12, 
and 15 (FPM 10, FPM 13, and WMB 55), the positions are from the 3.5 cm map 
of Tr 13. For source 13, the position is from the 6.2 cm map of Tr 14. 

^ Source numbers from IRS: Taylor & Storey 1984, B: Brown 1987, WMB: 
Wilking at al. 1997, FCW: Feigelson et al. 1998, CT: Choi & Tatematsu 2004, 
SMM: Nutter et al. 2005 and Groppi et al. 2007, FPM: Forbrich et al. 2006, and 
SMA/IRAC: Groppi et al. 2007. 



pected for the free-free emission from coUimated thermal radio 
jets or coronal winds (—0.1 < a < 1; Reynolds 1986; Anglada 
et al. 1998). However, there are a few exceptions. The spectral 
index of B 9 is negative in both Tr 1 and Tr 13/14 estimates, 
suggesting that the centimeter emission is from a nonthermal 
process at least partially. The spectral index of IRS 7B needs 
some discussion (see § 5.2). The spectral index of B 5 varies 
wildly, which is probably due to a combined effect of nonther- 
mal emission processes and flare activities (see § 5.4). 

It should be noted that, even though the spectral index of IRS 
5 is consistent with free-free emission, its circularly polarized 
flux certainly tells that it has a nonthermal emission component. 
There are other examples of nonthermal emitter with a slightly 
positive spectral index, such as T Tau Sb (Loinard et al. 2007). 
These sources suggest that the spectral index alone is not a very 
reliable indicator of emission process. On the other hand, in ra- 
dio surveys of star forming clouds, sources displaying negative 
spectral indices are usually considered background nonthermal 
sources (e.g., Rodriguez et al. 1999), but the existence of young 
stellar objects such as B 9 tells that the spectral index alone 
is not a good criterion for distinguishing young stellar objects 
from extragalactic objects. 

4.3. Long-Term Variability 

Figure 2 shows long-term (8-13 years) variabiUty of the 3.5 
cm flux density. Considering the month-scale variability, most 
sources do not show clear trends. IRS 7A seems to show a 



steady increase of flux density. More observations are needed 
to confirm this trend. If real, such long-term variability is likely 
due to changes in outflow activity rather than circumstellar mag- 
netic activities. 

4.4. Prope r Motion 

Proper motions of bright compact sources were measured by 
comparing the uniform- weight 3.5 cm maps of Tr 1 and Tr 13. 
Most sources did not show a detectable proper motion. R CrA 
is the only source showing a marginally detectable proper mo- 
tion: 0.23" ± 0.11" toward the southeast over the 8.1 yr time- 
baseline, or (20, -21) mas yr^'^, which is consistent with the 
proper motion of the R CrA association, (5.5, -27.0) mas yr~^ 
(Neuhauser et al. 2000). 

4.5. Polarimetry 

Circularly polarized 3.5 cm emission was detected toward 
three sources while none of the radio sources showed detectable 
polarization in the 6.2 cm band. IRS 5b showed detectable po- 
larization in most observing runs, but IRS 5a was never detected 
in the Stokes-F flux maps. Polarized emission from IRS 7A and 
B 5 was detected once each. Details of the polarimetry results 
will be presented separately in Paper II (Choi et al. 2008). 

4.6. Infrared Photometry 

Since flux measurements and images based on Spitzer data 
are reported previously for many of the sources in the R CrA re- 
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Fig. 1. — Finding cliarts of the radio continuum sources in tlie R CrA region, (a) Map of tlie A = 3.5 cm continuum. This 3.5 cm map was made by 
averaging all the images from tracks Tr 1-13. For each track, a natural-weight image was made. Each of these images was corrected for the primary 
beam response and convolved to an angular resolution of 3", and then all the 13 images were averaged. The contour levels are 1, 2, 4, 8, 16, 32, 
and 64 times 0.09 mjy beam"'^. (b) Natural- weight maps of the A = 3.5 and 6.2 cm continuum toward the IRS 7 region observed in 2005 February 
(Tr 13/14). The contour levels are 1, 2, 4, 8, 16, 32, and 64 times 0.06 mJy beam~^. Dashed contours are for negative levels. The rms noise is 0.02 
mJy beam^^. Shown in the bottom right corner are the synthesized beams. The straight line near the bottom left corner corresponds to 1000 AU at 
a distance of 170 pc. Plus signs: Submillimeter continuum sources SMM IC, SMM IB, and SMM lA, from north to south (Groppi et al. 2007). 
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Table 3 

Flux Densities of the Radio Sources in the R CrA Region 





F{62 cm) 


F(3.5 cm) 




Source Name 


Peak 


Total 


Peak 


Total 


a 


Trl: 1996 December 


1 IRS 2" 


0.15 ± 0.03 


0.25 ± 0.04 








2 B5 


0.46 ± 0.03 


1.06 ±0.11 


0.42 ± 0.04 


0.51 ±0.07 


-1.3 ±0.3 


3 IRS Sa'^ 






0.11 ±0.02 


0.12 ± 0.02 




4 IRS Sb"^ 






0.17 ±0.02 


0.25 ± 0.02 




3+4 IRS 5^= 


0.22 ± 0.02 


0.43 ± 0.03 




0.35 ± 0.07 


-0.4 ±0.4 


5 IRS 6 












6 IRS 1 


0.16 ±0.02 


0.16 ± 0.02 


0.22 ± 0.02 


0.26 ± 0.03 


0.9 ±0.3 


7 R CrA 


0.10 ±0.02 


0.11 ±0.04 


0.19 ±0.02 


0.19 ±0.02 


1.0 ±0.7 


8 FPM 10 






0.08 ± 0.02 


0.11 ±0.02 




9 B 9b'*'= 


0.39 ± 0.03 


0.32 ± 0.03 


0.15 ±0.02 


0.19 ±0.04 


-0.9 ±0.4 


10 B 9a<^ 


0.40 ± 0.03 


1.09 ± 0.03 


0.49 ± 0.02 


0.90 ± 0.04 


-0.35 ± 0.09 


9+10 B 9 




1.33 ± 0.03 




1.09 ± 0.04 


-0.36 ± 0.08 


11 IRS7A 


2.61 ± 0.02 


3.63 ± 0.07 


2.40 ± 0.02 


3.71 ± 0.07 


0.04 ± 0.05 


12 FPM 13 


0.10 ±0.02 


0.16 ±0.03 


0.08 ± 0.02 


0.14 ±0.02 


-0.2 ±0.4 


13 CHLT 13 












14 IRS 7B 


0.92 ± 0.02 


2.54 ±0.11 


0.73 ± 0.02 


1.16 ±0.03 


-1.41 ±0.09 


15 WMB 55 












16 T CrA 


0.13 ±0.02 


0.13 ±0.02 


0.10 ±0.02 


0.10 ±0.02 


-0.5 ±0.5 


Tr 13/14: 2005 February 


1 IRS 2'' 


0.17 ±0.03 


0.17 ±0.03 








2 B5 


1.06 ± 0.03 


1.21 ± 0.04 


1.78 ± 0.05 


2.18 ±0.14 


1.06 ±0.13 


3 IRS 5a'=''= 






0.16 ±0.03 


0.16 ±0.03 




4 IRS 5b'^''= 






0.41 ± 0.03 


0.41 ± 0.03 




3+4 IRSS'^ 


0.58 ± 0.02 


0.77 ± 0.03 




0.77 ± 0.05 


0.00 ±0.14 


5 IRS 6 












6 IRS 1 


0.27 ± 0.02 


0.30 ± 0.04 


0.46 ± 0.02 


0.50 ± 0.05 


0.9 ±0.3 


7 R CrA 


0.15 ±0.02 


0.15 ±0.02 


0.28 ± 0.02 


0.28 ± 0.02 


1.1 ±0.3 


8 FPM 10 


0.27 ± 0.02 


0.42 ± 0.03 


0.19 ±0.02 


0.30 ± 0.04 


-0.6 ±0.3 


9+10 B 9' 


1.27 ±0.02 


1.76 ± 0.07 


1.14 ±0.02 


1.56 ± 0.04 


-0.22 ± 0.09 


11 IRS7A 


5.08 ± 0.02 


6.24 ± 0.08 


4.22 ± 0.02 


6.09 ± 0.06 


-0.04 ± 0.03 


12 FPM 13 


0.31 ±0.02 


0.47 ± 0.04 


0.16 ±0.02 


0.52 ± 0.06 


0.2 ±0.3 


13 CHLT 13 


0.17 ±0.02 


0.17 ±0.02 








14 IRS 7B 


0.89 ± 0.02 


2.56 ± 0.12 


0.81 ± 0.02 


2.31 ± 0.15 


-0.19 ±0.14 


15 WMB 55 






0.15 ±0.02 


0.17 ±0.03 




16 T CrA 






0.13 ±0.02 


0.17 ±0.02 




Note. — Flux densities 


are corrected for the primary beam response, and the unit is mjy (beam For 


most sources, the flux densities are from the natural-weight maps. In most cases, total flux of each source was 


measured in a box, typically a circumscribed rectangle of la contour. 






^ Spectral index between 6.2 cm and 3.5 cm. In 1996 December, the 6.2 cm and the 3.5 


cm observations 


were made in a single track. In 2005 February, the 6.2 cm observations were made 9 days after the 3.5 cm 


observations, and the spectral index of rapidly varying sources should be interpreted with caution. 


IRS 2 is located outside the field of view of the 3.5 cm maps. 






° IRS 5a and 5b were not clearly separated in the 6.2 cm maps. 






The 6.2 cm fluxes of B 9a and 9b are from the uniform-weight maps. 






The 3.5 cm fluxes of these sources are from the uniform- 


■weight maps. 






^ B 9a and 9b were not clearly separated in either the 6.2 cm or the 3.5 cm maps. 
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Fig. 2. — Variations of 3.5 cm flux densities. Filled circles: Flux densities from this work. Open circles: Rux densities from previously published 
works (Suters et al. 1996; Feigelson et al. 1998; Forbrich et al. 2006, 2007; Miettinen et al. 2008). 



Table 4 

Infrared Flux Densities 



Source 


5(3.6 //m) 


5(4.5 //m) 


5(5.8 //m) 


5(8 ^lm) 


5(24 //m) 


5(70 ^lm) 


IRS 2 


2.57 


3.65 


6.95 


9.61 


174 


35.7 


IRS 6 


0.099 


0.162 


0.21 


0.23 


0.44 


4.7 



Note. — Fluxes are from the Spitzer data, and the imit is Jy. Photometric uncertainties are ~10%. 
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gion (Groppi et al. 2007; Forbrich & Preibisch 2007), we report 
the photometry of IRS 2 and IRS 6 only (Table 4). The IRAC 
color-color diagram for the infrared sources in the R CrA region 
is shown in Figure 3. 

For IRS 2 and IRS 6, the bolometric temperature and lumi- 
nosity were calculated using the Spitzer data along with the Two 
Micron All Sky Survey (2MASS) data. Although IRS 2 has a 
greater bolometric luminosity than IRS 6 by a factor of 25, they 
have a similar bolometric temperature of ^450 K, indicative of 
Class I sources (Chen et al. 1995). According to the [3.6]-[4.5] 
versus [5.8]-[8.0] color-color diagram (Allen et al. 2004; Lee et 
al. 2006), both IRS 2 and IRS 6 are also classified as embedded 
Class O/I sources (Fig. 3). 



DISCUSSION 



5.7. 



IRS 7 A Region 

The IRS 7A region is a very complicated site of star forma- 
tion. There are at least two young stellar objects in this small re- 
gion, IRS 7A and B 9 (SMA 2), and each of them have subcom- 
ponents. All these sources seem to be embedded in the dense 
core SMM ICQ (Nutter et al. 2005; Groppi et al. 2007). Fig- 
ures 4-6 show the centimeter continuum maps of this region. 
We will try to dissect this region by comparing the radio maps 
with the images in other wavelengths. 

^In this paper, we use the submillimeter source numbers given by Nutter et al. 
(2005). Choi & Tatematsu (2004) referred to SMM lA, IB, and IC as vdA 4, 
5, and 3, respectively, following the source numbers given by van den Ancker 
(1999). There is a difference in astrometry between Nutter et al. (2005) and 
Groppi et al. (2007), and we use the positions given by Groppi et al., if possible. 



RS 7E 



IRS 7A 



Class 



Class 0/1 



Class 



IRS 5N 
WMB 55 

- IRS 5 

IRS 6 
RS 1 

- IRS 2 

- R CrA 

- T CrA 

- IRS 4 



[5^8]-[8.0] 

Fig. 3. — Color-color diagram for the young stellar objects in the R 
CrA region from the Spitzer IRAC 3.6, 4.5, 5.8, and 8.0 /im data. For 
IRS 2 and IRS 6, the photometric data are listed in Table 4. For the 
other sources, the flux densities are from Groppi et al. (2007). The 
dotted boxes delineate the domains of Class O/I, II, and III sources 
(Allen et al. 2004; Lee et al. 2006). 
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Fig. 4. — Map of the 3.5 cm continuum toward the IRS 7A region from 
high-resolution data sets: Tr 1, 2, 3, 4, and 13. For each data set, a 
natural-weight image was made. Each of these images was corrected 
for the primary beam response and convolved to an angular resolution 
of 1.2" (shown in the bottom right corner), and then all the five im- 
ages were averaged. The contour levels are 1, 2, 4, 8, 16, 32, 64, and 
128 times 0.025 mjy beam^ ^ . Dashed contours are for negative levels. 
Squares: The 7 mm continuum sources CT 2, 3, 4, and 5, from north 
to south (Choi & Tatematsu 2004). Plus sign: The 1.1 mm continuum 
source SMA 2 (Groppi et al. 2007). 




5 0-5 -10 

RA offset (arcsec) 



5 0-5 -10 

RA offset (arcsec) 



Fig. 5. — Uniform-weight maps of the 6.2 cm continuum toward the 
IRS 7A region. The contour levels are 1, 2, 4, 8, 16, and 32 times 0.10 
mJy beam~^. Markers are the same as those in Fig. 4. Observation 
dates and track numbers are labeled. Shown in the bottom right corner 
are the synthesized beams (see Table 1). 
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Fig. 6. — Uniform-weight maps of the 3.5 cm continuum toward IRS 7A and B 9. The contour levels are 1, 2, 4, 8, 16, and 32 times 0.06 mjy 
beam~^. A position correction was applied to the Tr 2 and 3 maps (see § 3.2). Plus sign: The 1.1 mm continuum source SMA 2 (Groppi et al. 2007). 
Squares: The 7 mm continuum sources CT 2, 3, and 4 (Choi & Tatematsu 2004). Cross: The X-ray source CXOU J190155.3-365722 (Hamaguchi 
et al. 2005). Star symbol: The mid-IR source IRAC 5 (Groppi et al. 2007). 
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5.1.1. IRS 7A 

The Tr 1 map of Figure 6 has the highest resolution and 
shows the accurate position of the compact component of IRS 
7A. In addition to the compact component, the other maps of 
Figure 6 as well as the average map (Fig. 4) show the extended 
structure that seems to be elongated roughly in the north-south 
direction. Judging from the spectral index and the morphology, 
the centimeter emission of IRS 7A most likely traces a thermal 
radio jet. The outflow activity may have become enhanced in 
recent years as the 3.5 cm flux density in 2005 (Tr 13) is sig- 
nificantly higher than those in the 1996-1998 period (Fig. 2). 
The 3.5 cm map of Tr 13 (Fig. 6) also shows a newly appearing 
subcomponent at 0.85" southeast of the main peak, which may 
be a bright knot of recently ejected material. Since this subcom- 
ponent is not present in the images of earlier epoch, its proper 
motion is at least 130 mas yr^^ (^100kms~^ at 170 pc). 

The peak position of the 3.5 cm maps has not changed over 
the 8 years of monitoring and may indicate the location of the 
central object. As shown in Figure 6, the 3.5 cm source posi- 
tion agrees well with the positions of the sources detected in 
other wavelengths: the mid-lR source IRAC 5, the 7 mm source 
CT 4, and a Chandra X-ray source (Groppi et al. 2007; Choi 
& Tatematsu 2004; Hamaguchi et al. 2005). Despite all these 
detections over a wide range of wavelength, the classification 
(hence the evolutionary status) of IRS 7A is frustratingly un- 
certain because of the complicated environment. IRS 7A is cer- 
tainly not a T Tauri star The nondetection at 1 . 1 mm (Groppi 
et al. 2007) seems to rule out the possibility of IRS 7A being 
a Class source, and Hamaguchi et al. (2005) suggested that 
IRS 7A may be a Class I source seen at a large inclination angle 
based on the near-IR magnitude and the column density mea- 
surement from X-ray observations. The IRAC color-color dia- 
gram places IRS 7A near IRS 7B and above the other Class 1 
sources (Fig. 3), which suggests that IRS 7A may be either a 
Class O/I transitional object or a Class I object in the younger 
part of the class. 

Hamaguchi et al. (2005) could reproduce the XMM-Newton 
spectrum of the X-ray source Xw using a single-temperature 
plasma model, which could also reproduce the Chandra data. 
Comparing the 3.5 cm maps with the Chandra X-ray positions 
(Fig. 6), the X-ray source seems to be associated with the main 
radio peak itself but not with the subcomponent on the south- 
east side or the more extended component on the south side of 
the main peak. (The Chandra observations were made in 2000 
October and 2003 June.) This geometry suggests that the X-ray 
source may be closely related with the central star. The X-ray 
emission may come from either the protostellar magnetosphere 
or near the base of outflow. The detection of circularly polarized 
radio emission (Paper 11) strongly suggests the existence of an 
active magnetosphere. 

5.1.2. B 9 Complex 

The centimeter source B 9 was resolved into two sources, B 
9a and 9b, with a separation of 0.5" (Fig. 6). B 9a is slightly 
elongated in the east-west direction with a deconvolved size of 
~0.5", and B 9b is a point-hke source. Both of them show neg- 
ative spectral indices (-0.9 to -0.2, see Table 3), suggesting that 
the centimeter emission comes from a nonthermal process at 
least partially. In contrast, Miettinen et al. (2008) reported pos- 
itive (0.4) spectral index, which suggests that either the non- 
thermal component is highly variable or their flux measurement 
included a significant contribution from the extended thermal 
emission source. The proper motion of B 9ayb is uncertain be- 



cause they were not clearly separated in the Tr 13 image. 

The centimeter source B 9 and the 7 mm source CT 2 are 
located close to each other, but not exactly at the same position. 
It is unlikely that this position difference (1.2") is caused by an 
instrumental uncertainty because the position agreement of IRS 
7A and CT 4 is good. The peak position of the 1 . 1 mm source 
SMA 2 is 0.4" south of B 9a, but B 9b and CT 2 are also located 
within the beam size of the 1.1 mm map of Groppi et al. (2007). 
We will use the name 'B 9 complex' to refer to all these compact 
sources collectively. The spectral energy distribution (SED) of 
the B 9 complex (see Fig. 8 of Groppi et al. 2007) suggests that 
it may be harboring a Class source, but it is not clear what is 
the true nature of the B 9 complex. Considering the complex- 
ity of the centimeter-millimeter maps and the existence of two 
outflows (see below), we suggest that the B 9 complex may be 
a multiple-protostar system. A subarcsecond-resolution map of 
dust emission is needed to locate the central stellar object(s) and 
imderstand the situation better. 

Radio sources with significantly negative spectral index are 
rare in star forming regions but not unique. Examples in- 
clude outflow lobes of Serpens SMM 1 and IRAS 16547^247 
(Rodriguez et al. 1989; Garay et al. 2003). The suggested ra- 
diation mechanism is nonthermal synchrotron emission, but the 
existence of highly relativistic electrons is not easy to explain. 
It is possible to accelerate a small fraction of electrons in the 
shocked region of outflow (e.g., Crusius-Watzel 1990), but the 
density of the medium should be low enough to avoid being 
dominated by thermal radiation (Henriksen et al. 1991). It is 
also possible to produce energetic electrons in active magneto- 
spheres, as observed around more evolved young stellar objects 
(Massi et al. 2006, 2008; Loinard et al. 2007). Other types of 
nonthermal radiation (such as gyrosynchrotron radiation, elec- 
tron cyclotron maser, and plasma maser) cannot be ruled out 
(Giidel 2002), but they are less likely because circular polariza- 
tion was not detected toward B 9 (Paper 11). It remains to be 
seen if B 9a/b are indeed shocked regions of outflowing gas, 
and measurements of proper motion and hneiir polarization will 
be helpful. 

5.1.3. FPM 10 and FPM 13 

FPM 10 and 13 are extended sources, and their shape and 
peak position change depending on the epoch of observations. 
Their centimeter spectra are flat, suggesting optically thin free- 
free emission. They are highly elongated (Figs. 1, 4, and 5), 
and the major axes point to the B 9 complex. Therefore, they 
are most likely ionized gas of outflows driven by young stellar 
object(s) in the B 9 complex. However, they cannot be a single 
bipolar outflow system because their directions with respect to 
B 9 are nearly perpendicular to each other Forbrich & Preibisch 
(2007) found diffuse X-ray emission in a region just north of 
FPM 13, but the relation between the extended X-ray structure 
and the B 9 complex is not clear. 

FPM 10 is located about 3" west of the B 9 complex and 
elongated in the east-west direction. The size along the ma- 
jor axis is ~3". The location and the flow direction of FPM 
10 match well with the redshifted lobe of the bipolar outflow 
mapped in several molecular lines including HCO+ J = 1^ 
and CO J = 3 ^ 2 (Anderson et al. 1997; Groppi et al. 2004, 
2007). Therefore, FPM 10 and the molecular bipolar outflow 
seem to belong to a single outflow system. 

FPM 13 is located about 4" north of the B 9 complex and 
elongated in the north- south direction. The size along the major 
axis is ~3". There is no known molecular outflow that can be 
associated with FPM 13, but this is probably because the molec- 
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ular line maps do not have a sufficiently high angular resolution. 
Because of the flow direction, the driving source of FPM 13 may 
be an object in the B 9 complex but cannot be the same object 
that is driving FPM 10. It should also be noted that we caimot 
completely rule out the possibility that FPM 13 may be driven 
by IRS 7A. 

5.1.4. CT3 

CT 3 was nearly as bright as CT 2 in the 7 mm map (Choi 
& Tatematsu 2004) but not clearly detected in the centimeter 
maps. It is located in the "bridge" of radio emission structure 
between IRS 7A and B 9 (Fig. 4). This "bridge" may be tracing 
an outflow either from IRS 7A or from the B 9 complex (or a 
superposition of outflows from both). A marginally detectable 
amount of 3 .5 cm emission can be seen in the Tr 1 3 map (Fig. 1 b 
left panel). Choi & Tatematsu (2004) suggested that CT 3 may 
be a deeply embedded young stellar object, but the nondetection 
at 1 . 1 mm (Groppi et al. 2007) seems to rule out this possibility. 
A possible explanation is that CT 3 may be a knot of outflowing 
gas that was temporarily brightened at the time of the 7 mm 
observations. 

5.2. IRS 7B 

The IRS 7B region is less complicated than the IRS 7A re- 
gion but still needs a careful analysis to understand the struc- 
ture. Figures 7-9 show the centimeter continuum maps of the 
IRS 7B region. There is a single compact source surrounded by 
an extended structure elongated in the northeast-southwest di- 
rection. The spectral index was highly negative (-1.4) during 
Tr 1 but only mildly negative (-0.2) during Tr 13/14 (Table 3). 
This difference is due to the low 3.5 cm flux during Tr 1, which 
probably means that the extended structure was resolved out in 
the A-array observations. If only the compact source is consid- 
ered, the spectral index is ^0.1 during Tr 1 and ^0.7 during Tr 
13/14. Therefore, the centimeter flux from the central compact 
source of IRS 7B may be mostly from free-free emission. 

The peak position of the 3.5 cm maps did not show a de- 
tectable proper motion over the 8 years of monitoring and may 
indicate the location of the central object. As shown in Figure 9, 
the 3.5 cm peak position agrees with the positions of the mid-IR 
source IRAC 4, the 1.1 mm source SMA 1, and the Chandra X- 
ray source (Groppi et al. 2007; Hamaguchi et al. 2005). Groppi 
et al. (2007) suggested that IRS 7B is a Class 0/1 source based 
on the SED analysis. Therefore, IRS 7B may be a protostar em- 
bedded in the dense core SMM IB. 

The extended component of IRS 7B shown in Figure 7 has a 
length of ^9" and a position angle of ~20°. The bipolar mor- 
phology suggests that the extended component traces a bipo- 
lar outflow driven by IRS 7B. The northeastern lobe (FPM 15) 
is brighter than the opposite lobe. The outflow in our maps is 
probably the base of the large-scale (~2') bipolar jet in the 6 
cm maps of Brown (1987) and Miettinen et al. (2008). Large- 
scale maps of molecular lines from single-dish observations did 
not show any sign of molecular outflow parallel to the IRS 7B 
outflow. Again, this nondetection is probably caused by a confu- 
sion with the east- west outflow of B 9. The interferometric map 
in the HCO+ line (Fig. 6 of Groppi et al. 2007) seems to show 
a molecular outflow around SMA 1: blueshifted emission on 
the south of SMA 1 and redshifted emission on the northeast of 
SMA 1. The size and direction of the HCO"'" emission structure 
seem to agree with the outflow seen in our centimeter contin- 
uum maps. However, note that Groppi et al. (2007) interpreted 
the HCO+ emission structure as a part of a large spherical shell. 



u 

M 



I 



o 

I 



IRS7B 



CD 




3.5cm 




15 



10 



RA offset (arcsec) 



Fig. 7. — Map of the 3.5 cm continuum toward the IRS 7B region from 
high-resolution data sets: Tr 1, 2, 3, 4, and 13. See Fig. 4 for details. 
The contour levels are 1, 2, 4, 8, 16, and 32 times 0.025 mjy beam~^. 
Plus sign: The 1.1 mm continuum source SMA 1 (Groppi et al. 2007). 
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Fig. 8. — Uniform-weight maps of the 6.2 cm continuum toward the 
IRS 7B region. The contour levels are 1, 2, and 4 times 0.10 mJy 
beam~^. Plus sign: The 1.1 mm continuum source SMA 1 (Groppi 
et al. 2007). Star symbol: The mid-IR source IRAC 4 (Groppi et al. 
2007). 
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Fig. 9. — Uniform-weight maps of the 3.5 cm continuum toward IRS 7B. The contour levels are 1, 2, 4, 8, and 16 times 0.06 mjy beam ^. Plus 
sign: The 1.1 mm continuum source SMA 1 (Groppi et al. 2007). Cross: The X-ray source CXOU J190156.4-365728 (Hamaguchi et al. 2005). 
Star symbol: The mid-lR source IRAC 4 (Groppi et al. 2007). 



Similarly to IRS 7A, the Chandra X-ray source is associated 
with the main peak of the centimeter source, IRS 7B, but not 
with the secondary peaks in the outflow lobe such as FPM 15. 
This geometry suggests that the X-ray source is associated with 
the central protostar. Hamaguchi et al. (2005) found that the 
XMM-Newton spectra of Xe could be reproduced better using 
a two-temperature model. They concluded that the hard X-ray 
component comes from a mass accretion spot on the protostar 
while the soft component probably comes from shock-heated 
plasma related to a jet/outflow. The high plasma temperature 
and the large X-ray luminosity require a large amount of gas at 
a very high outflow velocity (Hamaguchi et al. 2005; Forbrich 
et al. 2006). A possible model reconciling the radio and the X- 
ray observations is that the hard and time-variable X-ray com- 
ponent may be emitted by circumstellar magnetic activity and 
that the soft and constant X-ray component and the centimeter 
emission may come from plasma heated by the collision of a 
steady jet with ambient gas (Hamaguchi et al. 2005; Forbrich & 
Preibisch 2007). Since the X-ray spectra from Chandra obser- 
vations were dominated by the soft component and the spatial 
resolution of XMM-Newton is limited (~6"), the exact position 
of the hard component is unknown. To understand the IRAS 7B 
system clearly, it is necessary to observe the X-ray source with 
a high resolution during an active phase of the hard component. 



5.2.1. CHLT13 

CHLT 13 is a 6.2 cm continuum source detected during Tr 
14 (Fig. \h). Since it was not detected in the 6.2 cm map of Tr 
1 and never detected in the 3.5 cm continuum, it must be a tran- 
sient source. Its location is near the dense core SMM 1 A. Nutter 
et al. (2005) suggested that the submiUimeter source SMM lA 
may be a prestellar object because it was not detected in other 
wavelengths. On the other hand, the 6 cm maps of Miettinen et 
al. (2008) show extended clumpy structure around this region, 
but none of their emission peaks coincides with CHLT 13 ex- 
actly. On the nature of CHLT 13, we suggest two possibilities. 
First, CHLT 13 may indicate a star forming activity in SMM 
lA. Second, CHLT 13 may be a temporarily brightened knot of 
the IRS 7B outflow. More speculatively, the southwestern out- 
flow of IRS 7B may be colliding with a dense clump in SMM 
1 A and producing transient radio emission at the shock front. 

5.3. IRS 5 

IRS 5 is a very interesting binary system. Figures 10-12 
show the centimeter continuum maps of the IRS 5 region. IRS 
5 was clearly resolved into two sources when the resolution is 
high enough (see the maps of Tr 1, Tr 2, and Tr 13 in Fig. 
12). Even when the resolution is not high enough, the 3.5 cm 
maps usually show an extended structure that can be best in- 
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Fig. 10. — Map of the 3.5 cm continuum toward the IRS 5 region from 
high-resolution data sets: Tr 1, 2, 3, 4, and 13. A natural-weight image 
was made from the Tr 1 data set, and a uniform- weight image was made 
from each of the other data sets. Each of these images was corrected 
for the primary beam response and convolved to an angular resolution 
of 0.9" (shown in the bottom right comer), and then all the five images 
were averaged. The contour levels are 1, 2, 4, 8, and 16 times 0.04 
mjy beam^^. Plus sign: The submillimeter continuum source SMM 
4 (Groppi et al. 2007). Open star symbols: The mid-IR sources IRAC 
9 (IRS 5N) and IRAC 10 (Groppi et al. 2007). Filled star symbols: 
Near-IR sources listed in the 2MASS Point Source Catalogue. 
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IRS 5 region. The contour levels are 1, 2, 4, and 8 times 0.07 mJy 
beam^^. Star symbols: The same as those in Fig. 10. Crosses: Sources 
detected in X-rays (Hamaguchi et al. 2008). 



terpreted as a combination of two point-like sources (see the 
maps of Tr 4 and Tr 6-12 in Fig. 12). In addition, comparison 
with the Stokes-1/ maps (Paper 11) clearly shows that the south- 
western source was detected in the Stokes-/ maps. Since the 3.5 
cm sources are similar to the near-lR sources (Chen & Graham 
1993; Nisini et al. 2005) in separation and position angle, we 
identify the southwestern and the northeastern sources as IRS 
5a and 5b, respectively. The Stokes-F maps show a point-like 
source at the position of IRS 5b while IRS 5a was undetected in 
the Stokes-F flux, and there is a good correlation between the 
Stokes-y and Stokes-/ fluxes of IRS 5b (see the discussion in 
Paper 11). 

IRS 5a is brighter in the near-lR images of Nisini et al. 
(2005), but IRS 5b is usually brighter in radio. The separation 
between them in the 3.5 cm Tr 1 image is 0.9" (160 AU at 170 
pc). They were not clearly separated in the 6.2 cm maps, but 
the Tr 1 map (Fig. 1 1 left panel) shows that both IRS 5a and 5b 
have detectable 6.2 cm fluxes. Since circularly polarized flux 
was detected from IRS 5, the centimeter flux clearly has a non- 
thermal component, and we will discuss this issue in Paper 11. 
Another Class 1 source in this region, IRS 5N (FPM 5), was not 
clearly detected though Forbrich et al. (2006) reported that it 
was detected when several data sets were combined. 

Table 5 lists the 3.5 cm flux densities of each binary com- 
ponent. For the maps of Tr 4-12, the flux was measured by fit- 
ting the image with a sum of two point-like sources (see the 
end note of Table 5 for details). The uncertainty of total flux 
from this fitting method was estimated by performing simple 
Monte Carlo simulations. Since we assumed that both sources 
are point-like, the uncertainty can be larger if this assumption is 
not valid. The light curves of IRS 5a and 5b are shown in Fig- 
ure 13. The flux uncertainty shown in Figure 13 is the statistical 
uncertainty listed in Table 5, but one should also consider the 
uncertainty in the absolute flux scale, which was not explicitly 
measured but usually about 1-2% (VLA Calibrator Manua0)- 
The radio light curves clearly show that both sources are highly 
variable. IRS 5b is usually brighter and more active than IRS 5a, 
but IRS 5a outshone its companion at least once (Tr 4) during 
the monitoring. This behavior clearly tells that the flux variabil- 
ity of each source should be analyzed separately, but previous 
authors considered the flux of the IRS 5 system as a whole, 
which may have caused some confusion. 

The radio positions of IRS 5a/b agree well with the position 
of X-ray sources (Hamaguchi et al. 2008). The separation be- 
tween the X-ray sources is ^^0.8", which is in good agreement 
with the radio value. There seems to be a difference in the rel- 
ative degree of flare activity between the two sources. IRS 5b 
is more active in radio (Fig. 13), but IRS 5a is more active in 
X-rays (Hamaguchi et al. 2008). To understand this difference, 
it is necessary to monitor the IRS 5 system in both radio and 
X-rays simultaneously. Forbrich et al. (2007) attempted such a 
study, but IRS 5 was not active in radio during their observations 
and their angular resolution was not high enough to resolve the 
system. 

The position and relative brightness of the near-lR counter- 
parts need some consideration. Nisini et al. (2005) reported that 
the separation between the two near-lR sources is ^0.6", but 
their Figure 1 shows that the separation is actually ^0.9", in 
agreement with the radio/X-ray value. The absolute position of 
the near-lR sources in the image of Nisini et al. (2005) has an 
offset of ~1" to the southeast relative to the radio/X-ray posi- 
tion. We presume that this offset was caused by an instrumental 

^ See |http://www. vla.Drao.edu/astro/calib/manual | 
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uncertainty of the near-IR observations. It is interesting to note 
that the 2MASS position agrees with the radio position of IRS 
5b while the mid-lR source IRAC 10 coincides with IRS 5a 
(Fig. 10). This difference between the 2MASS and the IRAC 
positions is not an astrometric problem because the positions 
of IRS 5N show a good agreement. Therefore, we suppose that 
the primary object in the near-IR brightness can be either one 
depending on the epoch of observations. That is, either one or 
both of them may be variable in near-IR. Alternatively, they may 
have different intrinsic colors. The relation between the variabil- 
ity in near-IR and in radio/X-rays has not been well understood 
so far (Forbrich et al. 2007). 

IRS 5a/b (and IRS 5N) are young stellar objects probably 



embedded in the dense core SMM 4 (Nutter et al. 2005; Groppi 
et al. 2007). IRS 5 as a whole is a Class 1 source (Fig. 3; Wilk- 
ing et al. 1997; Nutter et al. 2005; Hamaguchi et al. 2008). 
However, Nisini et al. (2005) suggested that IRS 5a may be a 
more evolved object, such as a T Tauri star, on the basis of the 
lack of significant accretion activity. It is not clear if IRS 5 has 
any outflow activity. Groppi et al. (2007) reported a CO out- 
flow in this region though the driving source could not be speci- 
fied clearly. Subarcsecond-resolution imaging in the millimeter 
or submillimeter band is necessary to understand the situation 
more clearly. 
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Table 5 

Flux Densities of the IRS 5 Binary Components at 3.5 cm 



Track 


IRS 5a 


IRS 5b 


Peak 


Total 


Peak 


Total 


Trl 


nil + f) 09 


19 + 00') 






Tr2 


1 + 03 


V/. u yJ.yJJ 


1 6^ ± n fi^ 


1 S? -1- 04 


Tr3 






W.J 1 31 V.VO 


U.Jo 31 U.UH- 


Tr4 








W.JJ zn u.ut- 


Tr5 










Tr6 




0.16 ± 0.03 


0.45 ± 0.03 


0.45 ± 0.04 


Tr7 




0.28 ± 0.04 


1.10 ±0.02 


1.05 ±0.04 


Tr8 




0.26 ± 0.05 


0.94 ± 0.02 


0.91 ± 0.03 


Tr9 




0.27 ± 0.04 


2.67 ± 0.02 


2.67 ± 0.04 


Tr 10 




0.41 ± 0.03 


1.77 ± 0.03 


1.83 ± 0.04 


Trll 




0.22 ± 0.04 


3.11 ±0.02 


3.14 ± 0.05 


Tr 12 




0.08 ± 0.04 


1.38 ±0.02 


1.43 ±0.04 


Tr 13 


0.16 ±0.03 


0.16 ±0.03 


0.41 ± 0.03 


0.41 ± 0.03 



Note. — Flux densities are corrected for the primary beam response, and the 
unit is mjy (beam~^). The Tr 1 fluxes are from the natural- weight map, and the 
others are from the uniform-weight maps. For Tr 1, 2, 3, and 13, IRS 5a and 5b 
are relatively well-separated in the image, and each total flux was measured in 
a box around each source. For Tr 4-12, the binary components were not clearly 
separated, and the total fluxes were measured by fitting each image using a sum 
of two point-like sources. In the fitting process, the position offset between the 
two sources was fixed using the value measured in the Tr 1 map. 



5.3.1. IRS 5b 

The flux variability of IRS 5b needs a detailed discussion. 
Though more extensive monitoring is necessary to understand 
the flare activities, it may be useful to consider two phases: qui- 
escent phase and active phase, depending on the total flux. We 
suggest a flux threshold of 0.8 mJy for IRS 5b, which is some- 
what arbitrary. The flux density in the quiescent phase (Tr 1, 3, 
4, 5, 6, and 13) is about 0.5 mJy on average. 

Examining the hght curve shown in Figure 13, two periods 
of active phase can be recognized: a short one in 1997 January 
(Tr 2) and a long one in 1998 September-October (Tr 7-12). 



3.5 - 

3 ^ 

2.5 J 

- 2 ^ 

1.5 ^ 

1 ^ 
0.5 

: • 

L 



1997 
JAN 



JUL 



1998 

AUG SEP OCT 
Trl 1 



Tr9: 



200^ 

FEB^ 



Tr2 



Tr4 



Tr7 



2450 



3000 3050 

Time (Julian Date — 24' 



3100 



5400 



Fig. 13. — Light curves of IRS 5a/b at 3.5 cm. Open circles: Flux den- 
sities of IRS 5a. Open triangles: Flux upper limits (3cr) for IRS 5a (Tr 
3/5). Filled circles: Flux densities of IRS 5b. 



Detailed information on the 1997 flare activity is not available, 
and the duration of the active phase must have been less than 
22 days. The 1998 flare activity can be discussed in more de- 
tail. The duration of the active phase was about 30 days, and 
the flux fluctuated wildly during this period. The average flux 
density during this active phase was higher than that of the qui- 
escent phase by a factor of about 4. The light curve shows three 
flux peaks (Tr 7, 9, and 11) with peak-to-peak intervals of about 
10 days, but the flux might have fluctuated more rapidly if the 
sampling interval had been shorter. 

The time scale of the flare activity of IRS 5b is much longer 
than that of the radio bursts of the Sun or flare stars, typically 
a few minutes or shorter. This month-scale flare activity of IRS 
5b implies that the cause of the flare is not a local magnetic 
activity on the surface of the (proto)star but probably a recon- 
nection event of global magnetic fields around the star or even 
magnetic fields connecting the star and the circumstellar disk. 
More details of the flare activity will be discussed in Paper II. 

5.3.2. IRS 5a 

IRS 5a is usually weak (~0.13 mJy) or undetectable, but 
showed at least two active periods during the monitoring. (Here, 
we set the quiescent-active phase threshold at 0.2 mJy.) In 1998 
June (Tr 4), IRS 5a went through an outburst with the flux 
higher than the quiescent level by a factor of ~15. It is not clear 
how long was the duration of this outburst. 

In 1998 September-October (Tr 7-1 1), the flux of IRS 5a in- 
creased mildly and stayed at a level about twice the quiescent 
level. This increase of flux is not an artifact of the image fitting 
process used in measuring the flux of each binary component. 
If the flux of IRS 5a had been at the quiescent level, it would 
have been barely noticeable in the images, but IRS 5a is clearly 
recognizable in the images during this active period. For exam- 
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Fig. 14. — Natural-weight maps of the 3.5 cm {left two panels) and 6.2 cm (right two panels) continuum toward B 5. The contour levels are 1, 2, 4, 
and 8 times 0.12 mjy beam~^. Plus sign: The 3.5 cm continuum peak position of the Tr 1 map. 



pie, compare the Tr 9 image with the Tr 12 image: the source 
structures in these two images are clearly different even though 
the beam sizes are almost identical, which clearly shows that the 
flux of IRS 5a was significantly higher during Tr 9 than during 
Trl2. 

5.3.3. Concurrent Activity in 1998 

What is really interesting about the flare activity in 1998 
September-October is that both objects in the binary system are 
active at the same time. If this concurrency is not a chance coin- 
cidence but has a real physical reason, there are three possibil- 
ities. First, the binary companions may be interacting directly 
through a physical process unknown so far. Second, what was 
observed may be a giant magnetic structure such as a circumbi- 
nary halo of ionized gas. Third, they may be responding to an 
external cause, such as an enhanced accretion of material from 
the envelope. 

The possibility of the interacting (proto)binary requires an 
underlying physical mechanism that is unknown so far. If the 
intense flare activity of IRS 5b was the cause of the flux en- 
hancement of IRS 5a, the light curves indicate that such an in- 
teraction should happen in a time scale of several days. The pro- 
jected separation between IRS 5a and 5b is ~0.9 hght-day. Then 
the order-of-magnitude velocity of the interaction agent must be 
about 10% of the speed of light or faster. Such a high velocity 
is not completely unreasonable because the circularly polarized 
radio emission of IRS 5b implies the existence of mildly rela- 
tivistic electrons. Intrabinary magnetic fields between stars are 
often suggested to explain the activity of stellar systems such as 
RS CVn binaries (Dulk 1985), and interacting helmet stream- 
ers were observed around a young binary system (Massi et al. 
2008). The size of the IRS 5 system is probably too large to have 
directly interacting stellar magnetic fields. However, large-scale 
magnetic fields involving circumstellar disk(s) and/or outflows 
might provide a new mechanism of interaction. 

Circumbinary halos and large-scale magnetic structures are 
often suggested to explain the radio activity of RS CVn bina- 
ries and related objects, and some halos were imaged using the 
VLBI technique (Mutel et al. 1985; Phillips et al. 1991). It is 
unknown if a circumbinary halo can be as large as the IRS 5 sys- 
tem. In principle, this possibility can be tested by imaging IRS 
5 with an angular resolution of ^0.2" during an active phase. 

In the third possibility, the responses of IRS 5a and 5b can 
be synchronized down to the crossing time of the accreting ma- 



terial. The velocity of accreting gas can be approximated by the 
free-faU velocity, vs « (2GM/r)V^. Assuming M = IMg 
and taking r to be the binary separation, we get vs 3 km 
s~^. Then the crossing time is ~200 yr, which is too long to ex- 
plain the concurrent activities of IRS 5a/b. Therefore, the third 
possibility is improbable unless there is a very fast triggering 
mechanism. 

5.4. B 5 

B 5 is a variable radio source (Fig. 2; Suters et al. 1996). 
Figure 14 shows the radio continuum maps of B 5. While the 
3.5 cm map of Tr 1 did not have enough sensitivity to see the 
source structure, the 6.2 cm map of Tr 1 shows that B 5 is 
extended along the northeast-southwest direction (RA. = 44°) 
with a deconvolved size (FWHM) of 1.3". In the 3.5 cm map 
of Tr 13, B 5 appears slightly extended along a similar direc- 
tion (RA. = 53°) with a deconvolved size of <0.6". The proper 
motion measurement of B 5 is unreUable because the source is 
extended. The peak position of the 6.2 cm source has a small 
offset (0.6" toward northeast) with respect to the position of the 
3.5 cm source. 

The spectral index was highly negative (-1.3) during Tr 1 
but positive (1.1) during Tr 13/14 (Table 3). This difference is 
due to the high 3.5 cm flux during Tr 13, which probably means 
that B 5 went through an outburst during Tr 13. This outburst 
was accompanied with circularly polarized emission, suggest- 
ing magnetospheric origin (Paper 11). The centimeter flux dur- 
ing Tr 1 was obviously from a nonthermal emission process. 
Such a large fluctuation of spectral index was also reported by 
Suters et al. (1996) based on monitoring observations in 1992. 
However, the spectral index was never highly negative in the 
1992 observations, and Suters et al. concluded that the fluctua- 
tion was caused by an optical depth effect of thermal free-free 
emission. Miettinen et al. (2008) reported a negative spectral in- 
dex, confirming the nonthermal nature of the emission from B 5 . 
The complicated variability of the centimeter emission from B 
5 may be caused by a combined effect of nonthermal emission 
process and frequent flare events. 

The nature of B 5 is an enigma. Suters et al. (1996) argued 
against an extragalactic source on the basis of the rapid vari- 
ability. The source size also supports this argument because the 
radio flux from such an extended source cannot vary rapidly if 
it is extragalactic. Examinations of archival Spitzer images did 
not produce any positive detection at the position of B 5. The 
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nondetection in the mid-IR or millimeter continuum suggests 
that it is not a deeply embedded object if located at the distance 
of the CrA cloud. Feigelson et al. (1998) suggested B 5 to be a 
brown dwarf based on a faint source in the near-lR images of 
Wilking et al. (1997), but this region is confusing in optical and 
near-IR images because of the nebulosity around IRS 2. 

5.5. IRS 1 

IRS 1 (HH 100 IR) is a point-like source, and its flux is rel- 
atively stable. Figure 15 shows the centimeter continuum maps 
toward the IRS 1 region. The spectral index is positive and con- 
sistent with optically thick thermal emission, which agrees with 
the conclusion of Suters et al. (1996). IRS 1 probably showed a 
proper motion of 21 ± 10 mas yr~^ toward the east, but more 
observations are needed to confirm this value. 

The 3.5 cm peak position agrees well with the position of the 
mid-IR source IRAC 7 (Groppi et al. 2007). IRS 1 is a variable 
X-ray source (Forbrich et al. 2006), which suggests that the low 
variability of the radio flux is probably due to the large optical 
depth in the centimeter band. IRS 1 is a Class I object embed- 
ded in SMM 3, and most of its bolometric luminosity is due to 
accretion (Fig. 3; Wilking et al. 1992; Nisini et al. 2005; Nutter 
et al. 2005; Groppi et al. 2007). It is the driving source of HH 
100 and probably other Herbig-Haro objects (Strom et al. 1974; 



Hartigan & Graham 1987; Wang et al. 2004). 

5.6. IRS 2 

IRS 2 is a point-like source, and Figure 16 shows the radio 
continuum maps. IRS 2 is located outside the field of view in 
our 3.5 cm observations, but is within the field of view of the 
3.5 cm observations of projects AC 481 and AM 596. The 3.5 
cm flux density varies between 0.25 and 0.64 mJy. The spec- 
tral index between 3.5 and 6.2 cm, estimated from the average 
flux densities, is 1.3 ± 0.7, but this value is highly uncertain 
because there was no observation made simultaneously in both 
wavelength bands. The positive spectral index implies that ther- 
mal free-free radiation may be the emission mechanism at the 
centimeter band. 

The 3.5 cm peak position agrees well with the near-IR source 
position (Fig. 16). IRS 2 is a strong and variable X-ray source 
(Forbrich et al. 2006). IRS 2 seems to be embedded in SMM 
5 and shows evidences of active accretion (Nisini et al. 2005; 
Nutter et al. 2005). The bolometric temperature and luminos- 
ity of IRS 2, calculated using the infrared fluxes from Spitzer 
archival data (Table 4), are 457 ± 17 K and 10.0 ± 0.1 Lq, 
respectively. This temperature confirms that IRS 2 is a Class I 
protostar as previously suggested by Wilking et al. (1992). The 
IRAC color-color diagram also places IRS 2 inside the Class O/I 
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domain (Fig. 3). 

5.7. RCrA 

R CrA is a point-like source, and Figure 17 shows the radio 
continuum maps. While R CrA is dominating this region in the 
optical and near-IR bands, it is relatively weak in radio. R CrA 
is detectable in all the data sets analyzed and shows some flux 
variability (Fig. 2). The spectral index of the centimeter con- 
tinuum is positive, suggesting that the radio flux comes from 
optically thick thermal emission. The proper motion of R CrA 
is consistent with that of the R CrA association (see § 4.4). 

The 3.5 cm peak position agrees well with the near-IR source 
position (Fig. 17). It is not clear how strong the submiUimeter 
emission from R CrA is (Groppi et al. 2007). R CrA is a Herbig 
Ae star of ^3.5 Mq, and it is also an interesting X-ray source 
with a plasma temperature of ~100 MK (Forbrich et al. 2006). 
To explain the X-ray emission, Forbrich et al. (2006) suggested 
that R CrA may be a binary system: an optically bright Herbig 
Ae star and an X-ray bright Class I source. Takami et al. (2003) 
also suggested that R CrA may be a binary system based on 
spectro-astrometric observations, but, in this case, both mem- 
bers must be visible in the Ha hne. 



5.8. T CrA 

Radio emission from T CrA has not been reported before, 
but it was detected in our maps as a weak point-like source. 
Careful examinations of archival data revealed that T CrA was 
detectable most of the time. Figure 18 shows the radio contin- 
uum maps, and Table 6 Usts the flux densities. The spectral in- 
dex is negative, but the uncertainty is too large to constrain the 
emission mechanism. 

The 3.5 cm peak position agrees well with the position of the 
mid-lR source IRAC 1 (Groppi et al. 2007). T CrA is a weak X- 
ray source (Forbrich et al. 2007) and was not detected in the 
submiUimeter continuum (Groppi et al. 2007). T CrA is a Her- 
big star of type FOe. Takami et al. (2003) suggested that T CrA 
may be a binary system. Wang et al. (2004) suggested that T 
CrA may be the driving source of the Herbig-Haro object HH 
733. 

5.9. WMB55 

Centimeter continuum emission from WMB 55 has not been 
reported before, but it was detected in our 3.5 cm map as a weak 
point-Uke source. Examinations of archival data revealed that 
the 3.5 cm flux of WMB 55 was detectable in a few observing 
runs, but it was not detected in the 6.2 cm continuum. Table 6 
lists the flux densities, and Figure 19 shows the 3.5 cm contin- 
uum maps. The nondetection in the 6.2 cm continuum implies 
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Table 6 

Flux Densities of T CrA and WMB 55 at 3.5 cm 





TCrA 


WMB 55 


Track 


Peak Total 


Peak Total 



Trl .. 
Tr2 . . 
Tr3 . . 
Tr4 . . 
Tr5 .. 
Tr6.. 
Tr7 . . 
Tr8 . . 
Tr9 . . 
Trio. 
Trll . 
Tr 12. 
Tr 13 . 



0.10 ±0.02 
0.11 ±0.02 
0.11 ±0.02 
0.10 ± 0.02 

0.19 ± 0.02 
0.11 ±0.03 
0.11 ±0.03 
0.16 ±0.03 
0.08 ± 0.02 
0.16 ± 0.03 



0.10 ±0.02 
0.21 ± 0.03 
0.11 ±0.02 
0.10 ± 0.02 

0.20 ± 0.02 
0.13 ±0.03 
0.13 ±0.03 
0.16 ±0.03 
0.08 ± 0.02 
0.16 ± 0.03 



0.13 ± 0.02 0.17 ±0.02 



0.13 



± 0.02 0.13 ± 0.02 



0.15 ± 0.03 0.15 ± 0.03 
0.18 ±0.03 0.20 ±0.04 
0.15 ±0.02 0.17 ±0.03 



Note. — Flux densities are corrected for the primary beam response, and the 
unit is mjy (beam~^). They are from the natural-weight maps. 
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Fig. 19. — Natural-weight maps of the 3.5 cm continuum toward WMB 55. The contour levels are 1 and 2 times 0.06 mJy beam ^. Plus sign: The 
submillimeter continuum source SMM 2 (Groppi et al. 2007). Star symbol: The mid-IR source IRAC 3 (Groppi et al. 2007). 



that the spectral index is positive, but the uncertainty is too large 
to constrain the emission mechanism. 

The 3.5 cm peak position agrees well with the position of 
the mid-IR source IRAC 3 (Groppi et al. 2007). WMB 55 was 
not detected in the X-ray band. WMB 55 is a Class I source 
embedded in SMM 2 (Fig. 3; Chini et al. 2003; Nutter et al. 
2005). Wang et al. (2004) suggested that WMB 55 is a probable 
driving source of the Herbig-Haro objects HH 733 and/or HH 
734. HH 104 A/B are also located in the vicinity of WMB 55, 
but it is not clear if they are related with WMB 55. 

5.70. IRS 6 

IRS 6 was not detected in our maps, though Forbrich et al. 
(2006) reported that it was detected when several data sets were 
combined. Surprisingly, examinations of archival data revealed 
that IRS 6 went through an outburst on 1998 October 13 (Tr 12) 
while it was undetectable or only marginally detectable in the 
other epochs. Miettinen et al. (2008) reported another outburst 
detected in 1998 January. 

Figure 20 shows the 3.5 cm continuum map at the epoch 
of the outburst. The source was slightly elongated along the 



northwest-southeast direction (RA. w -41 °) with a deconvolved 
size of 0.7", and the total flux density was 0.44 ± 0.05 mJy. 
Since IRS 6 was undetected during Tr 1 1, the onset time-scale 
was less than 3 days. It is not clear how long was the duration of 
this outburst because Tr 12 was the last observing track of the 
1998 observing campaign. 

The peak position at the outburst agrees well, within 0.2", 
with the X-ray source position, which suggests that the same 
source is responsible for both the radio and the X-ray emission. 
(The X-ray source position was determined using the images 
discussed in Hamaguchi et al. 2008.) Since IRS 6 is a near-IR 
binary system (Nisini et al. 2005), the radio/X-ray source may 
be one of the binary members. Comparing the near-lR source 
positions with the radio source position, after applying the as- 
trometric correction determined using IRS 5a/b (see § 5.3), the 
radio source can be identified as IRS 6a. Interestingly, the mid- 
IR source IRAC 8 is located halfway between IRS 6a and 6b, 
suggesting that the mid-IR fluxes of the binary members are 
comparable. X-ray emission was never detected at the position 
of IRS 6b. In contrast to the IRS 5 system, only one member of 
the IRS 6 system seems to show detectable magnetic activities. 
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Fig. 20. — Natural-weight map of the 3.5 cm continuum toward IRS 6. 
The contour levels are 1 , 2, and 4 times 0.06 mjy beam" ^ . Cross: X-ray 
source (Hamaguchi et al. 2008). Squares: Near-IR sources in Fig. 1 of 
Nisini et al. (2005). A position correction was applied using the offset 
determined with IRS 5a/b (see text). The binary separation is ^1.2". 
Open star symbol: The mid-IR source IRAC 8 (Groppi et al. 2007). 
Filled star symbol: Near-IR source listed in the 2MASS Point Source 
Catalogue. Note that this 2MASS source is not IRS 6 but a nebulosity. 
IRS 6 itself is listed in the Point Source Reject Table and coincides 
with IRAC 8 within ~0.2". 



IRS 6a showed little sign of accretion activity (Nisini et al. 

2005) and was suggested to be a T Tauri star (Forbrich et al. 

2006) . However, the spectral class of IRS 6 was unclear be- 
cause a detailed SED was not available. The 850 fim continuum 
map of Groppi et al. (2007) does show an unnumbered local 
peak at the position of IRS 6, which indicates that IRS 6 may 
be an embedded object. The bolometric temperature and lumi- 
nosity of IRS 6, calculated using the infrared fluxes from the 
Spitzer archival data (Table 4), are 424 ± 12 K and ~0.4 Lq, 
respectively. This temperature suggests that IRS 6 is a Class I 
protostar. The IRAC color-color diagram also places IRS 6 in 
the Class O/I domain (Fig. 3). Wang et al. (2004) suggested that 
IRS 6 may be driving a giant Herbig-Haro flow: HH 730, HH 
860, HH 104 C/D, and HH 99. 

6. SUMMARY 

We observed the R CrA region using the VLA in the 3.5 
and 6.2 cm continuum to image the young stellar objects and 
study the star forming activities in this region. Archival VLA 
data from recent observations made in extended configurations 
were also analyzed. A total of 16 sources was detected. The 
source properties were studied by investigating source structure, 
flux variability, spectral slope, and proper motion and also by 
comparing with the sources detected in shorter wavelength. The 
main results are summarized as follows: 

1 . The centimeter emission of IRS 7A traces the thermal ra- 
dio jet that seems to be flowing in the north- south direction, and 
the outflow activity has become enhanced in recent years. A re- 



cently ejected outflow knot was found to be moving at a speed of 
1 00 km s ^ ^ or higher. The main peak position of the radio emis- 
sion agrees with the mid-IR and X-ray positions, which suggests 
that the X-ray emission comes from the central (proto)star or a 
very vicinity of it. 

2. The Class source SMA 2 is associated with the centime- 
ter continuum sources B 9a/b and the 7 mm continuum source 
CT 2. While the detailed structure of the B 9 complex needs fur- 
ther studies, it may be a multiple-protostar system that is driv- 
ing at least two outflows (FPM 10 and FPM 13). The negative 
spectral slope of B 9a/b indicates that the radiation process is 
nonthermal. 

3. The centimeter continuum images of IRS 7B show a com- 
bination of a compact source and an extended structure. The 
compact source corresponds to the Class O/I source SMA 1 . The 
X-ray source position also agrees weU with the radio peak po- 
sition, and IRS 7B (SMA 1) is one of the youngest known ob- 
jects that were detected in X-rays, which suggests that energetic 
magnetic activities start very early in the protostellar stage. The 
extended structure is probably a bipolar outflow driven by IRS 
7B. 

4. The Class I binary system IRS 5 was resolved in the 3 .5 cm 
continuum maps. The binary separation is 0.9" (160 AU). Both 
members of the system show radio flare activities and X-ray 
emission. IRS 5 is one of the youngest known binary systems 
that is showing magnetic activities in both members. 

5. The radio light curve of IRS 5b shows at least two flare 
activities, and the duration of active phase is about a month. 
Such a long time-scale suggests that the flare activity involves 
large-scale magnetic fields. 

6. During flie 1 998 September-October flare event of IRS 5b, 
the radio flux of IRS 5a also showed a significant enhancement. 
This concurrent activity suggests a few possibilities: IRS 5a and 
5b are interacting direcdy, they are surround by a circumbinary 
halo, or they are responding to an external trigger. IRS 5 may 
be the youngest known example of interacting binaries. 

7. The variable radio source B 5 was found to be an extended 
(~1") object. The spectral slope also varies significantly, sug- 
gesting that both thermal and nonthermal radiation processes 
are working. The rapid variability and the source size rule out an 
extragalactic source. The nature of B 5, however, still remains 
unclear. 

8. The Class I sources IRS 1 and IRS 2 show positive spectral 
slopes. The centimeter emission of these sources is likely com- 
ing from optically thick thermal free-free radiation. The flux 
density of IRS 1 is relatively stable. 

9. The Herbig Ae star R CrA is a relatively weak radio 
source, and the spectral slope is consistent with thermal emis- 
sion. The marginally detected proper motion, 28 mas yr^^ to- 
ward southeast, is consistent with that of the R CrA association. 

10. Radio detections of the Herbig star T CrA and the Class 
I source WMB 55 are reported here for the first time. The emis- 
sion mechanism of these sources is not clear. 

1 1 . The near-lR binary IRS 6 is usually weak or undetectable 
in the centimeter continuum, but a radio flare event was de- 
tected in 1998 October. The radio source positionaUy coincides 
with the X-ray source and was tentatively identified as IRS 6a. 
Only one member of the IRS 6 binary system seems to show 
detectable magnetic activities. Infrared data suggests that IRS 6 
is a Class I source. 
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